Detection of focal source and arrhythmogenic substrate from body surface potentials to guide atrial fibrillation ablation

Focal sources (FS) are believed to be important triggers and a perpetuation mechanism for paroxysmal atrial fibrillation (AF). Detecting FS and determining AF sustainability in atrial tissue can help guide ablation targeting. We hypothesized that sustained rotors during FS-driven episodes indicate an arrhythmogenic substrate for sustained AF, and that non-invasive electrical recordings, like electrocardiograms (ECGs) or body surface potential maps (BSPMs), could be used to detect FS and AF sustainability. Computer simulations were performed on five bi-atrial geometries. FS were induced by pacing at cycle lengths of 120–270 ms from 32 atrial sites and four pulmonary veins. Self-sustained reentrant activities were also initiated around the same 32 atrial sites with inexcitable cores of radii of 0, 0.5 and 1 cm. FS fired for two seconds and then AF inducibility was tested by whether activation was sustained for another second. ECGs and BSPMs were simulated. Equivalent atrial sources were extracted using second-order blind source separation, and their cycle length, periodicity and contribution, were used as features for random forest classifiers. Longer rotor duration during FS-driven episodes indicates higher AF inducibility (area under ROC curve = 0.83). Our method had accuracy of 90.6±1.0% and 90.6±0.6% in detecting FS presence, and 93.1±0.6% and 94.2±1.2% in identifying AF sustainability, and 80.0±6.6% and 61.0±5.2% in determining the atrium of the focal site, from BSPMs and ECGs of five atria. The detection of FS presence and AF sustainability were insensitive to vest placement (±9.6%). On pre-operative BSPMs of 52 paroxysmal AF patients, patients classified with initiator-type FS on a single atrium resulted in improved two-to-three-year AF-free likelihoods (p-value < 0.01, logrank tests). Detection of FS and arrhythmogenic substrate can be performed from ECGs and BSPMs, enabling non-invasive mapping towards mechanism-targeted AF treatment, and malignant ectopic beat detection with likely AF progression

Challenges associated with interpreting mechanisms of AF

Determining optimal treatment strategies for complex arrhythmogenesis in AF is confounded by the lack of consensus regarding the mechanisms causing AF. Studies report different mechanisms for AF, ranging from hierarchical drivers to anarchical multiple activation wavelets. Differences in the assessment of AF mechanisms are likely due to AF being recorded across diverse models using different investigational tools, spatial scales and clinical populations. The authors review different AF mechanisms, including anatomical and functional re-entry, hierarchical drivers and anarchical multiple wavelets. They then describe different cardiac mapping techniques and analysis tools, including activation mapping, phase mapping and fibrosis identification. They explain and review different data challenges, including differences between recording devices in spatial and temporal resolutions, spatial coverage and recording surface, and report clinical outcomes using different data modalities. They suggest future research directions for investigating the mechanisms underlying human AF

Subject-Specific Ablation of Pathologic Conduction Patterns Beyond the Pulmonary Veins: A Personalised Modelling Approach

Improving patient outcomes with ablation of non-paroxysmal AF (PsAF) has proved challenging using a population-based treatment approach due to large interindividual variability in the underlying electroanatomical substrate. Ablation of pathologic conduction patterns outside of pulmonary vein isolation (PVI) has recently shown encouraging results in PsAF patients returning for their first or second retreatment (76% freedom from AF recorded in the RECOVER AF trial). However, the optimal targets and best sequence of ablation lesions are still unknown, and testing different sequences, types, and methods of ablation cannot be performed clinically on a single patient or patient cohort. Considering the predictive potential of computational modelling, a small exploratory subset of patients (N=4) enrolled in the ongoing DISCOVER trial was used to create patient-specific models of left atrial electrophysiology. The subject-specific models displayed a high correlation between simulated targets and clinical targets. AF complexity was highest in all patients prior to therapy. PVI caused a marginal decrease in complexity across the cohort whereas PVI+PCP showed an extensive decrease in the AF complexity across the patients and resulted in AF termination in all patients

Constructing a Human Atrial Fibre Atlas

Atrial anisotropy affects electrical propagation patterns, anchor locations of atrial reentrant drivers, and atrial mechanics. However, patient-specific atrial fibre fields and anisotropy measurements are not currently available, and consequently assigning fibre fields to atrial models is challenging. We aimed to construct an atrial fibre atlas from a high-resolution DTMRI dataset that optimally reproduces electrophysiology simulation predictions corresponding to patient-specific fibre fields, and to develop a methodology for automatically assigning fibres to patient-specific anatomies. We extended an atrial coordinate system to map the pulmonary veins, vena cava and appendages to standardised positions in the coordinate system corresponding to the average location across the anatomies. We then expressed each fibre field in this atrial coordinate system and calculated an average fibre field. To assess the effects of fibre field on patient-specific modelling predictions, we calculated paced activation time maps and electrical driver locations during AF. In total, 756 activation time maps were calculated (7 anatomies with 9 fibre maps and 2 pacing locations, for the endocardial, epicardial and bilayer surface models of the LA and RA). Patient-specific fibre fields had a relatively small effect on average paced activation maps (range of mean local activation time difference for LA fields: 2.67-3.60 ms, and for RA fields: 2.29-3.44 ms), but had a larger effect on maximum LAT differences (range for LA 12.7-16.6%; range for RA 11.9-15.0%). A total of 126 phase singularity density maps were calculated (7 anatomies with 9 fibre maps for the LA and RA bilayer models). The fibre field corresponding to anatomy 1 had the highest median PS density map correlation coefficient for LA bilayer simulations (0.44 compared to the other correlations, ranging from 0.14 to 0.39), while the average fibre field had the highest correlation for the RA bilayer simulations (0.61 compared to the other correlations, ranging from 0.37 to 0.56). For sinus rhythm simulations, average activation time is robust to fibre field direction; however, maximum differences can still be significant. Patient specific fibres are more important for arrhythmia simulations, particularly in the left atrium. We propose using the fibre field corresponding to DTMRI dataset 1 for LA simulations, and the average fibre field for RA simulations as these optimally predicted arrhythmia properties

Left atrial enhancement correlates with myocardial conduction velocity in patients with persistent atrial fibrillation

Background: Conduction velocity (CV) heterogeneity and myocardial fibrosis both promote re-entry, but the relationship between fibrosis as determined by left atrial (LA) late-gadolinium enhanced cardiac magnetic resonance imaging (LGE-CMRI) and CV remains uncertain. Objective: Although average CV has been shown to correlate with regional LGE-CMRI in patients with persistent AF, we test the hypothesis that a localized relationship exists to underpin LGE-CMRI as a minimally invasive tool to map myocardial conduction properties for risk stratification and treatment guidance. Method: 3D LA electroanatomic maps during LA pacing were acquired from eight patients with persistent AF following electrical cardioversion. Local CVs were computed using triads of concurrently acquired electrograms and were co-registered to allow correlation with LA wall intensities obtained from LGE-CMRI, quantified using normalized intensity (NI) and image intensity ratio (IIR). Association was evaluated using multilevel linear regression. Results: An association between CV and LGE-CMRI intensity was observed at scales comparable to the size of a mapping electrode: −0.11 m/s per unit increase in NI (P < 0.001) and −0.96 m/s per unit increase in IIR (P < 0.001). The magnitude of this change decreased with larger measurement area. Reproducibility of the association was observed with NI, but not with IIR. Conclusion: At clinically relevant spatial scales, comparable to area of a mapping catheter electrode, LGE-CMRI correlates with CV. Measurement scale is important in accurately quantifying the association of CV and LGE-CMRI intensity. Importantly, NI, but not IIR, accounts for changes in the dynamic range of CMRI and enables quantitative reproducibility of the association

PIEMAP: Personalized Inverse Eikonal Model from cardiac Electro-Anatomical Maps

Electroanatomical mapping, a keystone diagnostic tool in cardiac electrophysiology studies, can provide high-density maps of the local electric properties of the tissue. It is therefore tempting to use such data to better individualize current patient-specific models of the heart through a data assimilation procedure and to extract potentially insightful information such as conduction properties. Parameter identification for state-of-the-art cardiac models is however a challenging task. In this work, we introduce a novel inverse problem for inferring the anisotropic structure of the conductivity tensor, that is fiber orientation and conduction velocity along and across fibers, of an eikonal model for cardiac activation. The proposed method, named PIEMAP, performed robustly with synthetic data and showed promising results with clinical data. These results suggest that PIEMAP could be a useful supplement in future clinical workflows of personalized therapies.Comment: 12 pages, 4 figures, 1 tabl